On the long-time behavior of spin echo and its relation to free induction decay 
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It is predicted that (i) spin echoes have two kinds of generic long-time decays: either simple 
exponential, or a superposition of a monotonic and an oscillatory exponential decays; and (ii) the 
long-time behavior of spin echo and the long-time behavior of the corresponding homogeneous free 
induction decay are characterized by the same time constants. This prediction extends to various 
echo problems both within and beyond nuclear magnetic resonance. Experimental confirmation of 
this prediction would also support the notion of the eigenvalues of time evolution operators in large 
quantum systems. 
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Calculations of spin echo (SE) envelope observed by 
nuclear magnetic resonance (NMR) in solids are, in gen- 
eral, very difficult to do and, in practice, mostly limited 
to simple sudden-jump models (e.g.,p|), which cannot 
capture complex non-Markovian dynamics induced by 
spin-spin interactions. The goal of this work is to ap- 
proach the SE problem from a new direction by making 
the following two propositions: 

1. The generic long-time behavior of SE envelope 
S(2t) obtained from the sequence [ir/2-pulse — time t 

— TT-pulse — time r — detection] has form either 

S(2t) = Ce- 2 ^ T (1) 

or 

5(2r) = e-^id + C 2 cos(2ujt + <f> s )], (2) 

where C, C\, C2, 7, u) and <f>s are some constants. 

2. The constants 7 and u> describing the long-time SE 
decay are the same as the long-time constants of the cor- 
responding "homogeneous" free induction decay (FID) 
defined by the pulse sequence [ir/2-pulse — time t — de- 
tection] and exhibiting generic long-time behavior |2j of 
the form either 



F(t) = Ae~^ 
counterpart of Eq. , or 

F(t) = Ae~ 7 *cos(wt + (j> F ). 



(3) 



(4) 



- counterpart of Eq.J2J). Here A and <pF are some ad- 
ditional constants. (The meaning of the "homogeneous" 
FID will be specified below.) 

An important consequence of the above statement is 
that the long-time constants (7 and uS) of SE envelope 
can be obtained from a much simpler FID calculation. At 
present, a number of approximations exist (see, e.g., 0,0 
and references therein), which can be used to calculate 
the FID long-time constants with reasonable accuracy. 

In the following, I shall first consider the case of spin- 
spin interaction only, and then generalize the treatment 



to the presence of spin-lattice interaction and to other 
echo problems both within and beyond NMR. 

The simplest setting, which exposes the generic dif- 
ficulties of the SE problem, involves two spin species / 
and S characterized by different Larmor frequencies. The 
relevant Hamiltonian is H = Jiinh + Ti-hom, where Tiinh 
contains the interaction of spins with a slowly varying 
in space external field; and Tthom includes the spin-spin 
interaction and, sometimes, the interaction with micro- 
scopically inhomogencous external fields. 

The SE technique cancels the effect of Hinh, so that 
the SE envelope measured on /-spins, is given by: 



S(2t) 



Tr 



iH h 



iHh 



where Ii X is the ^-component of the ith J-spin, and X(tt) 
is the operator of the 7r-pulse acting on /-spins. 

The process, which, according to Eq.(J5}, generates a 
point of the SE envelope can be expressed as follows: 



M(t) = 



F(t) 



0<t<T 
T <t, 



F(T,t) 

where F(t) is the homogeneous FID, i.e 



F(t) = Trie 



lH ho 



-in ha 



(6) 



(7) 



and F(t, t) is the "after-pulse" response of the system in 
the absence of inhomogeneous broadening (see Fig.^). 

Derivation of Eas. (|ll2[) involves two steps: step 1 — to 
show that the generic long-time behavior of FID is given 
by Eas.Q3l4|): and, step 2 — to show that Eqs.© and gj 
for FID imply Eqs.(UJ and © for SE. 

Step 1 has actually been made in Refs.0, 0] for 
the problems involving only one kind of equivalent 
spins. These references contain detailed arguments and 
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FIG. 1: An example of the process, which generates a single 
point (open circle) of the spin echo envelope [see Eq.@]. 



also summarize quite substantial experimental^ and 
numerical [(Ht} evidence that the long-time behavior (l.'il II) 
is generic for isolated classical and quantum spin systems. 
That argumentation applies directly to the present case, 
when Jihom is invariant under lattice translations. [The 
random placement of spins will be discussed later.] 

I now proceed to step 2 for the monotonic case 11 loll . 
The FID having the long-time behavior can be writ- 
ten as: 

„, . f non-universal : < t < To /cA 
F{t >~\Ae-^ : T <i, (8) 

where To is the characteristic time beyond which the 
difference between the actual FID and its exponential 
asymptotics becomes negligible. 

One can now observe that an exponential decay ap- 
pearing in an equilibrium correlation function [such as 
F(t)] likely indicates that, in the corresponding linear 
response non-equilibrium problem[2j, the leading correc- 
tion to the equilibrium density matrix of any subsystem 
of spins evolves as 



£>o(x)e 



-ft 



(9) 



where x is the set of variables describing the density ma- 
trix, and (?o(x) is some function. An analogous statement 
applies to the probability distribution of classical spins. 

In the case of classical hyperbolic chaotic systems, 
Eq.@ has a rigorous justification as it constitutes a 
manifestation of a Pollicott-Ruelle resonance 0, 0. Be- 
yond this very limited class of systems, that equation 
amounts to a very plausible hypothesis. In Refs.@,|j], the 
case has been made for extending the notion of Policott- 
Ruelle resonances to large isolated systems of classical 
and quantum spins. The numerical evidence for quan- 
tum analogs of Policott-Ruelle resonances has also been 
reported for kicked quantum top|lO| and kicked chain of 
spins 1/2 [H[l3. 

In the regime describable by Eq. © , the effect of iden- 
tical manipulation of the system (7r-pulse, in this case) 



should have self-similar effect on the density matrix up 
to the prefactor exp(— -ft). Since the same exponent also 
appears in Eq.JSJ , the after-pulse response for t > r 
can be factorized as 



F( T ,t) = F( T )F 1 (t-T), 



(10) 



where F i (t — t) is some function 

One can further stipulate that e -7 ' dominates the 
long-time behavior of F(t), not because of the specific 
initial conditions, but simply because (i) 7 is the eigen- 
value of a given time evolution operator; and (ii) among 
the eigenvalues selected by the symmetry of the initial 
conditions^], 7 has the smallest real part. Therefore, 
the density matrix created after the 7r-pulse should even- 
tually (after some delay time n ~ To) approach the same 
asymptotics, which implies that 



F 1 (t- T ) = 



non-universal 

B e -7(t-r) 



T < t < T- 

T + n <t. 



Tl 



(11) 



Thus, when r > max[T ,ri], Eas. H8llUlll|) yield 

S(2t) ee F(t, 2t) = F(t)F 1 (t) = ABe~ 2 ~< T , (12) 

which is equivalent to Eq.QJ. 

Now I repeat step 2 starting from the oscillatory FID 
In this case, 



where 



, , J non-universal : < t < To 
W= I /(*)+/*(*) : r <t, 



f(t) = ae-^ +l ^\ 



(13) 



(14) 



and a is a complex number. 

In the long-time regime, the density matrix should be 
dominated by two complex conjugate self-similar terms: 



Po(x)e 



-(7+iw)t 



(?o( x )e 



- (7— iui)t 



(15) 



where go (x) is now a complex- valued function. After the 
7r-pulse, 



F(r,t) = f(T)h(t -t) + f*(T)f*(t - r), 



(16) 



where 



fl(t-T) 



non-universal 



5 ie -(7+M(*-r) 



32e 



-(7-iw)(t-r) 



T < t < T + Tl 
T + Tl < t. 



(17) 

Here, 61 and 62 are complex numbers, which do not 
need to be the complex conjugates of each other, be- 
cause go( x ) is not real. Thus, when t > max[To,Ti], 
Ens.m 311411 (ill 711 give 

S(2t) = /(r)A(r) + /*(r)A*(r) = \a\e~^ (18) 
x [|&!|cos(2u;t - (j) a - (j) bl ) + |& 2 |cos(</> a + ^2)] , 
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where cf> a , 4>bi and <pb2 are the phases of a, b\ and b 2 , 
respectively. Equation l(TB|) is equivalent to Eq.J2J). 
Several comments and generalizations are now in order: 

1. When a non-Hamiltonian (spin-lattice) component 
is added to the above problem, it should only strengthen 
the reliability of Eas. (|3l4[) . because the Markovian pro- 
cesses explicitly imply time-irreversibility and generically 
lead to the exponentially decaying FIDs. The rest of the 
proof then remains unchanged. 

2. The exact study of echo in the presence of a sin- 
gle Markovian fluctuator^| fully supports the long-time 
asymptotics (111 11 . 

3. As generic, I consider the case, when ho ~ Wo, 
where ho is the typical value of local fields in Tihom, and 
Wo their characteristic fluctuation rate. In this case, To ~ 
Ti ~ W ~ . The "generic" deviation from the generic case 
can be described as a situation, when the delay times To 
and T\ are anomalously long. Such a situation may occur, 
when ho 3> Wo , or when the local fields or the fluctuation 
rates are broadly distributed. 

4. Of particular concern can be the case, when the lat- 
tice sites are randomly partitioned between spins / and S 
representing, e.g., two different nuclear isotopes. In this 
case, different /-spins are not equivalent to each other. 
Extending the analysis of Ref. 2], one can argue that, for 
each sub-ensemble of equivalent /-spins the long-time be- 
havior of FID is still given by Eas. (|3l4| ). and, therefore, 
the resulting signal F(t) is the superposition of a contin- 
uous distribution of exponents 1|3I4[1 . For nuclear spins 
in solids, the distribution of parameter 7 should be lim- 
ited from below, because every spin experiences a finite 
fluctuating local field. Therefore, as a matter of princi- 
ple, the long-time asymptotics 1MI 11 will, eventually, be 
reached. However, in the worst case scenario, the dis- 
tribution of 7 can be very broad, and, therefore, it can 
take impractically long time before the asymptotics (|3I4|) 
becomes pronounced. 

In many NMR settings, two factors, work against the 
above worst case scenario. Namely: (i) The long-range 
nature of magnetic dipolar interaction between nuclei im- 
plies that each /-spin interacts with a sufficiently large 
number of iS-spins and other /-spins. Therefore, differ- 
ent /-spins, which, in a strict sense, are not equivalent, 
experience statistically equivalent environments, which 
implies that the distributions of 7 and to are narrowly 
peaked, (ii) /-spins belonging to different sub-ensembles 
are still coupled to each other, which should make spin 
dynamics more "homogeneous" . 

It should be noted, however, that, the averaging of the 
oscillatory SE (j2Jl even over a narrow distribution of u> 
will, eventually, suppress the beats leaving one with the 
monotonic asymptotics 

5. Beyond NMR, the language of free induction de- 
cays and echoes is frequently applied to the situations, 
when there exists a distribution of local fields and their 
fluctuation rates. The arguments presented in this work 
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FIG. 2: Distribution of the local field fluctuation rates P(W) 
consistent with the present theory. 



apply to the cases, when the distribution of the fluctua- 
tion rates satisfy the condition shown in Fig. [21 Namely, 
the local fields can be divided in two groups: very slowly 
fluctuating — the effect of these fields is cancelled by 
the SE pulse sequence; and those fluctuating with rates 
of the order of 1 /to or higher — they control both the 
[homogeneous] FID and the echo envelope. 

6. Hu and Hartmann|l4l| have considered a mathemat- 
ical model of infinitely dilute subsystem of /-spins sur- 
rounded by a large density of randomly placed S'-spins in- 
teracting with /-spins via magnetic dipolar interaction. 
They have found that the long-time behavior of SE is 
given by e CT , where c is a constant. The difference 
between this asymptotics and the one given by Eq.Q 
is due to the fact that the above mathematical model 
allowed for a non-zero probability that spins S are ab- 
sent in any finite volume around spin /. Although the 
above probability decreases exponentially with the vol- 
ume, the corresponding fraction of /-spins produces an 
exponentially slow contribution to the SE decay. The 
competition between these two exponential factors then 
leads to a non-universal decay law. 

7. The experimental evidence for the monotonic echo 
asymptotics Q is easily available — see, e.g., Refs.fiH 
|l6(. The oscillatory SE decays are, presumably, more 
rare. Nevertheless, they have also been observed — see, 
e.g., Refs.[l7llili|. The data from Ref. 01 are reproduced 
in Fig. together with the fit of the form 

8. Everything in the above treatment of conventional 
SE applies to the situations, when the inhomogeneous 
broadening is absent, and the second pulse is not n but 
has a different "angle". This, in particular, includes the 
so-called solid echoes in one-component spin systems (but 
excludes the special pulse sequences designed to re- 
verse the spin dynamics [lSl|2fJ.) 

9. One can thus perform a conclusive experimental test 
of the present theory by studying solid echoes in CaF2. 
This well-characterized material exhibits the oscillatory 
FID asymptotics 5] . Such an experiment can verify not 
only the result of the theory (i.e. the echo envelope) but 
also the underlying picture, namely it can measure the 
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FIG. 3: Circles — experimental SE envelope for 95 Pt in Pt- 
Au alloy from Fig. 7 of Ref.[l7|. line — fit of the form 0. 
Inset: non- logarithmic plot of the same data. 




time, ms 

FIG. 4: Experimental data extracted from the oscilloscope 
photographs of 129 Xe transient signals presented in Fig.l of 
Ref . |21| : squares — FID; diamonds, triangles and stars are 
the solid echoes with the pulse sequences (in respective order) : 
90° — t — 90V , 90° - r - 180° o and 90° - r - 90° o . Solid 
lines are the exponential fits with the same time constant. 



process M(t) and check whether the long-time behavior 
of F(t, t) is characterized by the same time constants as 
F(t) independently of the type of the pulse applied at 

t = T. 

10. The latter kind of test can, in fact, be performed on 
the old results reported for the solid mixture of 129 Xe — 
131 Xe 0| . The experimental traces of 129 Xe FID and 
three different solid echoes are reproduced in Fig. 0] As 
seen in this figure, the long-time behavior of the FID and 
that of two echoes can be well fit by the same exponent. 
The third echo, however, does not approach clear expo- 
nential shape within the observational time window. A 
conclusive test of the theory would, obviously, require 
more accurate measurements. 

11. Echo studies |2^| have indicated that some kind 
of slow fluctuators of unknown microscopic origin are re- 
sponsible for decoherence in solid-state g-bits. It is possi- 
ble that, in each particular g-bit, a small number of such 
fluctuators control the decoherence process. If true, the 
fluctuation rates of those few fluctuators can be widely 
separated from each other and, therefore, divided into 
"slow" and "fast" as shown in Fig. [5] In such a case, the 
present theory is applicable and can be used to analyze 
the echo envelopes. In this context, I would like to point 
out that the long-time behavior @ is consistent with the 
beats in the echo envelope observed in Ref.|22|. 

In conclusion, the analysis presented in this work sug- 
gests that, when the long-time behavior of SE is observed 
to be of the form Q or (J2J), the time constants of this 
decay can be extracted from the corresponding homoge- 
neous FID calculation. The experimental confirmation 
of the present theory (in particular, the solid echo exper- 
iment in CaF2) would also give a strong support to the 
notion of the eigenvalues of time evolution operators in 
quantum systems. The author is grateful to CP. Slichter, 
J. Haase and M. Mehring for discussing this work. 
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